Characterization of Fluctuations of Impedance and Scattering Matrices in Wave 

Chaotic Scattering 
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In wave chaotic scattering, statistical fluctuations of the scattering matrix S and the impedance 
matrix Z depend both on universal properties and on nonuniversal details of how the scatterer is cou- 
pled to external channels. This paper considers the impedance and scattering variance ratios, VR Z 
and VR S , where VR Z = Var[Zij]/ ' {Var[Z ti \Var[Z ]] \} 1/2 , VR S = Var [Sij]/{Var[Su] Var [Sjj}} 1/2 , 
and Var[.] denotes variance. VR Z is shown to be a universal function of distributed losses within 
the scatterer. That is, VR Z is independent of nonuniversal coupling details. This contrasts with 
VR S for which universality applies only in the large loss limit. Explicit results are given for VR Z 
for time reversal symmetric and broken time reversal symmetric systems. Experimental tests of the 
theory are presented using data taken from scattering measurements on a chaotic microwave cavity. 
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I. INTRODUCTION 

The general problem of externally generated time har- 
monic waves linearly interacting with a structure of lim- 
ited spatial extent is basic in many fields of science 
U HI • In recent years much work has been done elu- 
cidating the consequences for the scattering of waves in 
cases in which, in the geometric optics approximation, 
the ray orbits within the structure are chaotic. Examples 
include optical 3, acoustic 0, microwave 00 EE! an d 
electronic cavities 0, |n| . In the case of complex or ir- 
regularly shaped enclosures that are large compared with 
a wavelength, small changes in the frequency and the 
configuration give rise to large changes in the scattering 
characteristics. This feature motivates treatments that 
are statistical in nature. In this regard random matrix 
theory El E3 has proven useful in predicting universal 
aspect of chaotic wave scattering problems in the cases of 
both time reversal symmetric systems (corresponding to 
matrix statistics of the Gaussian Orthogonal Ensemble, 
GOE) and time reversal symmetry broken systems (cor- 
responding to matrix statistics of the Gaussian Unitary 
Ensemble, GUE). 

Scattering problems can be characterized by the scat- 
tering matrix S which relates outgoing scattered wave 
amplitudes b to incoming waves a, via b = Sa. An al- 
ternative formulation is in terms of the impedance ma- 
trix Z . To illustrate the impedance description, consider 
an electromagnetic wave scattering problem in which N 
transmission lines labeled i = 1, 2, • • • , N of characteris- 
tic impedance Zqi are connected to a cavity. Let Vj and 
Ii represent the voltage and current on transmission line 
i as measured at a suitable reference plane. Then the in- 
cident wave dj and the reflected wave 6, may be expressed 



as a t = (V, + Z 0i Ii)/Z^ 2 , h = (V % - Z 0i li) / Z'J 2 . The 
impedance matrix Z relates the vector voltage to the vec- 
tor current, via V = ZI, and Z and S are related by Z = 



Zq / (1 - + S)Z i ', where 1 is the TV-dimensional 

identity matrix, and Zq — diag(Zoi, Zq2, • • • , Zqn). The 
impedance formulation is identical to the so called "R- 
matrix" , a formulation introduced by Wigner and Eisen- 
bud in nuclear-reaction theory in 1947, and further de- 
veloped in Refs. Hi Elm Ell- 

Statistical variations of the elements of Z and S due to 
small random variations in the scattering are of great in- 
terest. These statistics have two fundamental influences, 
(i) universal aspects described by random matrix theory, 
and (ii) nonuniversal aspects dependent upon the details 
of the coupling of input channels (e.g., transmission lines) 
to the scatterer. Our main result concerns the quantity, 



VR Z = 



Var[Zi. 



y/Var[Zu]Var[Zj 



(1) 



where Var[A], the variance of the complex scalar A, is 
defined as the sum of Uar[i?e^4] and Var[ImA\. Our 
result is of the form 



VR Z 




for GOE, 
for GUE, 



(2) 



*Also at Department of Electrical and Computer Engineering. 



where A is a parameter characterizing the losses within 
the scatterer. For example, in the case of an electromag- 
netic cavity, A = uj/(2QAuj), where to is the frequency of 
the incoming signal, Aui is the average spacing between 
cavity resonant frequencies near o>, and Q is the qual- 
ity factor of the cavity (Q = oo if there are no internal 
losses). The remarkable aspect of @ is that (A) de- 
pends only on the loss parameter and not on the nonuni- 
versal properties of the coupling to the cavity. Thus VR Z 
is a universal function of the loss A. The results for F\ 
and Fi (to be derived subsequently) are shown in Fig. ^ 
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FIG. 1: Vi? z versus the loss parameter A, as specified in 
Eq. and Eg. iTSI. 



II. IMPEDANCE VARIANCE RATIO 

In this section we will obtain the universal functions 
Ft o (A ) in Eq. (J2J. We adopt a formulation (e.g., see Refs. 
(22.1231) that incorporates the nonuniversal effects of the 
specific coupling geometry of input-output channels to 
the scatterer, combined with the random matrix theory 
for the universal aspects of the chaotic wave behavior 
within the scatterer. (In what follows, we use terminol- 
ogy appropriate to microwave experiments.) Beginning 
with the case of zero loss (A = 0), we have that, in the 
GOE case, the impedance matrix Z is described by the 
following statistical model [22^. 



J 



E A ( fc «) 



R^ 2 (k n )w n w^ Rr' z (k n ) 



1/2, 



(7) 



For A > 1, 



For A < 1, 



VR, = 



VR Z = 



for GOE, 
for GUE. 



for GOE, 
for GUE. 



(3) 



(4) 



A ratio similar to 
scattering matrix 5 1 , 

VR S 



can also be considered for the 



VarlSi 



(5) 



^Var[Su}Var[S 3 j} 

In contrast with ©, VR S in general depends on both 
the coupling to the cavity and on the loss parameter A. 
However, in the special case of high loss, A ^> 1, VR S 
becomes universal, 



VR. = 



for GOE, 
for GUE. 



(6) 



That is, VR S = VR Z for A ^> 1. Based on their electro- 
magnetics experiments, Fiachetti and Michelsen |18| have 
recently conjectured the universality of © in the GOE 
case. More generally, © follows from classic results of 
Hauser and Feshbach describing fluctuations in the cross 
section of inelastic neutron scattering [Tsj . and this re- 
sult has been obtained by Friedman and Mello [23 using 
the concept of maximization of information entropy, and 
by Agassi et al. \2l\ using a random-matrix model. The 
important point is that a universal result for VR S [i.e., 
Eq. ijBJ] applies only for A> 1, while the universal result 
for VR Z , Eq. J2J and Fig. ^ is for arbitrary A. 

In what follows we derive these results and test them by 
comparison with data obtained from scattering measure- 
ments on a chaotic microwave cavity. Section II derives 
the results for impedance variance ration, VR Z . Section 

III considers the scattering variance ratio, VR S . Section 

IV presents our experimental tests of the theory. 



Here k is the wave number corresponding to the incoming 
frequency, w n is a vector whose elements are real, inde- 
pendent, zero mean, Gaussian random variables of unit 
variance. The eigenvalues are randomly chosen with 
statistics appropriate to GOE. This is done by generat- 
ing ensembles of eigenvalues of random matrices and then 
scaling them such that the mean spacing between adja- 
cent eigenvalues near fc„ is A(fc n ). Here A(fc) is given by 
the Weyl formula for the mean eigenvalue spacing. For 
example, for the case of a 2D cavity A(fc) is indepen- 
dent of fc, A = Att/ A where A is the area of the cavity. 
The system-dependent part of the coupling is character- 
ized by the corresponding radiation impedance matrix 
Z r = R r + jX r . Note also that, since Z r depends only 
on local geometry near the entrances to the scattering re- 
gion, its frequency dependence is much slower than that 
of Z which depends on the geometry of the scattering re- 
gion assumed to be characterized by much longer length 
scales. The radiation impedance is the impedance seen 
at the reference plane of the input channel when waves 
that leave the channel propagate outward and are not 
reflected back to their input by the distant walls of the 
cavity. Thus the radiation impedance is the impedance 
seen when the distant walls are removed to infinity, or (as 
can be done in an experimental measurement |9() when 
the distant walls are lined with absorber. Hence the ra- 
diation impedance depends only on the local structure in 
the vicinity of the port coupling to the cavity and not on 
the shape or chaotic properties of the cavity. In the case 
of ports that are far apart, e.g., of the order of the cavity 
size, the off-diagonal elements of Z r are small and will be 
neglected. Thus we will take Z r to be a diagonal matrix 
with elements Z r i = R rl + jX r i. In the GUE case, the 
elements of w n in Q are complex with real and imagi- 
nary parts each individually having Gaussian statistics, 
w n m i s replaced by (where f denotes the con- 
jugate transpose), and the statistics of are now those 
appropriate to GUE. 

Equation Q is derived in |23- It results from a for- 
mal series expansion of the solution for Z in terms of 
eigenfunctions of the closed cavity. The local structure 



3 



of the eigenfunctions is then assumed to satisfy the so- 
called " random plane wave hypothesis" ; namely that in 
the semi-classical limit is large compared to the typ- 
ical dimensions of the structure), if ray trajectories in the 
closed cavity are chaotic, then eigenfunction n is statisti- 
cally similar to a random isotropic superposition of plane 
waves of wave number k with magnitude \k\ — k n . This 
gives Eq. 2%. We then further supplement Eq. 
with the known distribution function of k 2 determined 
from random matrix theory and the local mean spacing 
Afc 2 between adjacent eigenvalues. 

The universality of the result for VR Z (Fig. 1) can 
be shown from as follows. The impedance matrix Z 
defined in will have a mean and a fluctuating com- 
ponent. As shown in refs [22ll23l |. the value of the mean 
is determined by all terms in the sum and thus depends 
on the slow /c-dependence of R r (k) and A(fc), while the 
fluctuating component of Z is determined by terms in 
the sum in Eq. for which k 2 ~ k 2 . Thus if R r {k n ) 
and A(fc„) are approximately constant over a range of 
A; 2 values corresponding to many resonances, then the 
fluctuating part of Z will be universal after appropriate 
normalization. If we define Z to be the fluctuating part 
of Z , we have 



R, 



-1/2 



{k)ZR^' 2 (k) = - ] -Y. 



,t 



c, 



(7') 



where s n are the eigenvalues of a GOE random matrix, 
normalized to have unit mean spacing between adjacent 
eigenvalues near s n ~ 0. Thus ( is a universal nor- 
malized impedance matrix. If Z r is diagonal, we have 
Cy = ^-ri^ij- Hence the ratios VR Z defined by 



ri rj l J ' 

is the same for Z and £, and is, therefore, universal. 
(We emphasize, however, that this conclusion relies on 
Z r being diagonal.) 

The effects of distributed loss, such as losses due to 
conducting walls or a lossy dielectric that fills the cavity, 
can be simply incorporated in Eq. 0. Since modal fluc- 
tuations in losses are small when the modes are chaotic 
and the wavelength is short, we can construct a com- 
plex cavity impedance accounting for distributed loss by 
simply replacing k 2 in Eq. by fc 2 (l — jQ _1 ), where 
Q is the cavity quality factor. (In terms of the normal- 
ized impedance matrix £ defined in (7') we replace the 
denominator s„ by s„ — jk 2 /[QA(fc)].) For example, 



Z„ 



n=l 



R„A n w 2 



-E 



R, 



(i-jQ- 1 )- 

&nW 2 n k 2 /Q 



Ri 



+ 3 



+ jXu 

R ri A n w 2 n {k 2 n 



pie, the expected value of Xa is, 



M 



E\Xii\ = lira - , 



71=1 



< k\i) 



dw in f{w vn )w 2 „ 
RriA n (k n k 



dk 2 



■ dk 2 j 



(k 2 



{k 2 /Qf 



(9) 

where f(wi n ) is the probability distribution function 
(pdf) of Wi n and Pj is the joint pdf of the eigenvalues. 
Integrating over all kj, j ^ n, we express E[Xa] as an 
integral over the pdf of k 2 , Pi(k 2 ) = 1/(A„M), we con- 



sider the M — > oo limit and use < 
random variable w„ ■ 



1 for the Gaussian 



E[X«] = / dk 



R n (k n )(k 2 -k 2 )/7r 
Vn {k 2 -k 2 ) 2 + (k 2 /Q) 2 



= X ri {k). (10) 



The second equality in (|1(J|) relating E[Xn] to the radi- 
ation reactance requires Q > 1 and is analogous to the 
Kramers-Kronig relation. 

The second moment of Xa can be determined in a 
similar way by integrating over all j except j = t, s and 
using the joint distribution function P2(k 2 ,k 2 ) = [1 — 



{\k 2 - k 2 \/A)}/(MA) 2 , where 



\kf 



k 2 1 /A) is known 
Using the fact that 



from Random Matrix theory 
g goes to zero at large argument and assuming that the 
radiation resistance R r i{k n ) and the average spacing A„ 
vary slowly over the damping width k 2 /Q, we obtain 



E[Xl\ = 



{E[X U }} 2 + 



(w 2 f 



k 2 /Q 
k 2 \/A)(k 2 ~k 2 )(k 2 



R 2 

2 



7T 
k 2 



dkfdk 2 



k 2 ) 2 + {k 2 /Q) 2 }[{k 2 - k 2 ) 2 + (k 2 /Q) 2 } 
Combining Eq. (jl 1|> and Eq. (|1C)|> . we obtain 



R 2 
A 



1 



(11) 



(12) 



4 + (x/A) 2 

where A = k 2 /(QA). A similar moment evaluation can 
be carried out for Ra, as specified in Eq. 0, which 
yields the same expression as Eq. Ijl2|) for Var[Ru]. For 
GOE (the case we are now considering) we have that 
0, g(s) = P(s) - [/ S d(s')f(s') - l/2](df/ds), where 
f(s) = [(sin7rs)/(7rs)]. 

In order to obtain the variance ratio, we also apply 
the previous process to the off diagonal term Z^j, i =/= j, 



k 2 ) 



which, based on Eq. 0, is given by 



kl) 2 + (k 2 /Q) 2> 
(8) 



(Henceforth, we employ the notation A n = A(k n ).) Cal- 
culation of the moments of the impedance is facilitated 
by the fact that the eigenvalues and eigenfunctions in the 
chaotic cavities are statistically independent. For exam- 



Zij — 



+ 3 



R r iR r j A n Wi n Wj n k 2 j Q 
(k 2 /Q) 2 + (k 2 - k 2 n ) 2 

Awi n W ■7 7 (. 



(13) 



'■J- 



(k 2 /Q) 2 + {k 2 - k 2 ) 2 



Since Wi n and Wj n are independent, the first moments of 
Xij and Rij are both zero, and the variance is equal to 



4 



the second moment, 



Var[X k 



M 



R ri R rj Al(wf n )(w 2 jn ) 

[{kl - k 2 ) 2 + {k 2 /QW 



M—>oc 7T Z 

R r i Rfj 1 
A 27r' 



(14) 

The same result is obtained for Var[Rij]. Combining 
Eq. UU) with Eq. (H2J), we have Eq. with 



^=F l( A) = [3-2^^ (x) _±_ r 



\ (15) 



A similar calculation in the GUE case is facilitated by 
the simpler form of the function g{x) which is now given 
by g(x) = sm 2 (nx)/(TTx) 2 . We obtain 



VR Z 



F 2 (X) = [2-2 



Sill 7TX , 



(x/xy 



= [i 



1 - e 



-4?rA 



4ttA 



(16) 

We note that subsequent to a previous announcement 
of our work [24[ , the two-frequency correlation functions 
for the elements of the impedance and the scattering ma- 
trix have recently been calculated by Savin, Fyodorov 
and Sommers |25| . and are consistent with the preceding 
in the limit of zero frequency separation. 



III. SCATTERING VARIANCE RATIO 

We now consider the scattering matrix in the high loss 
limit, A » 1. For simplicity, we consider the case of 
two channels connecting to the scatterer, N = 2, and 
Z and S are 2x2 matrices. We note that a chaotic 
scattering process can be divided into a direct process 
and a delayed process, which leads to a separation of the 
mean part (equal to Z r ) and the fluctuating part SZ of 
the cavity impedance, Z = Z r + SZ . The fluctuating 
part SZ decreases as loss increases. Thus in the high 
loss limit, SZ <C Z r , which implies ^12,^21 <C Zn 1 Z22- 
(Recall, the mean parts of the off diagonal components 

]/2 

are zero.) We may now form S using S = Z (Z — 



Z )(Z + Z a 



Since the off diagonal terms of Z 



are small, the diagonal components of S are dominated 
by the diagonal components of Z. We then find for Su, 



Su = 



Z\l — Z01 _ (Z r i — Z01) + SZ11 
Z\\ + Zqi (Z r i + Zqi) + SZn 



= Sri + [" 



2Z> 



01 



(17) 



;}8Zii, 



l (Z rl + Z 01 ) 21 

where S r i = {Z r \ — Zoi)/{Z r i + Zqi), and Zqi is the 
characteristics impedance of channel 1. Thus, we obtain 



Var[S! 



2Z i 



{Zri + Z m f 



l Var[Zi 



(18) 



In addition, we can express S12 in the high damping limit 
as 



S12 — 



2Zi2 v / ZqiZ i 



(12 



2Zi2 V ZaiZt 



02 



(Zu + Zqi){Z22 + Zm) {Z r \ + Zq\){Z t 2 + Zm) ' 

(19) 

which leads to 



Var[S 12 ] 



2V ZqiZq 



{Z r l + Zoi){Z r 2 + Z02) 



2 Var[Z 12 }, (20) 



and similarly for VarfS^i]. Combining Eq. (|18fl and 
Eq. 12()p . we recover Eq. © and we note that this re- 
sult is independent of the coupling (i.e., independent of 

Z r ). 

To illustrate the influence of coupling on VR S at finite 
loss parameter A, we consider the impedance matrix in 
the GOE case using the model normalized impedance 
C used in Ref. 22], Z = rI' 2 ^ 2 + jX r , where C is 
given by (y = -{j/n) Y,n=i(w t nW jn )/{k 2 - k 2 - jX), 
fe„ = k 2 1 A, and k 2 is set to be M/2, such that mean of 
£ is zero. Realizations of ( are produced numerically by 
generating Gaussian random variables Wi n and spectra 
kr\ from the eigenvalues of random matrices. We express 
a model scattering matrix S as 

s = U'^+n, + i)- 1 (7, 1/2 C7y 2 +.n, - 1), (21) 

where 7 r = Z^ 1 R r and j x = Z^' 1 X r , When j r is the 
identity matrix and "f x is zero, we reach the so-called 
perfect coupling condition, which means that the scat- 
tering is determined by the delayed process and the di- 
rect process is absent. We now consider an example 
in which the two-port couplings are the same so that 
Ir.x — diag{^ rxi ^ rx )i where % yX is a scalar. Figures 
121a) and (b) show results for the variation of VR S with 
the coupling parameters 7 r and j x , for a high loss case 
(A = 5) and for the lossless case (A = 0). In Fig. Eta), 
we fix j x to be zero, and vary 7 r , while in Fig. E^b), j r 
is fixed to be 1 and "f x is varied. Compared to the high 
damping case, VR S in the lossless case has a much larger 
deviation from the constant 1/2. Note that VR S is 1/2 in 
the perfect-coupling case (i.e., 7 r = 1, j x = 0), no mat- 
ter whether the cavity is highly lossy or lossless. This is 
related to the concept of "weak localization" reviewed in 

m 

In the case of an iV-port we can think of the above two- 
port consideration of VR S as applying to the iV-port con- 
verted to a two port by opening channels 3, 4, • ■ • , N; i.e., 
the incoming waves 03, 04, • • • aj\r are identically zero (for 
a microwave cavity with transmission line inputs, this 
corresponds to terminating transmission lines 3, 4, • • • , N 
with their characteristic impedances, Zq^, • • • , Zqn). 
Thus ports 3, 4, • • • , N effectively add to the loss due 
to the energy flux leaving through them. If the ports 
3, 4, • • • , N are assumed to act like distributed loss, they 
can be taken into account by increasing the loss pa- 
rameter A. [This increased loss enhances the validity of 
Eq. ©.] 
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FIG. 2: (a) VR S versus 7 r for 7^ = in the lossless case A = 
and in a high loss case A = 5. (b)VR s versus y x for j r = 1. 



IV. EXPERIMENTAL TESTS 

We provide experimental results testing the theoretical 
predictions for the statistical fluctuations in the variance 
of the <S* and Z elements, in the limit of large damp- 
ing. The experiments are done in an air-filled, quarter 
bow-tie shaped cavity which acts as a two-dimensional 
resonator below 19.05 GHz (Fig. 3(a)) This cavity 
has previously been used for the successful study of the 
eigenvalue spacing statistics 28] , eigenfunction statis- 
tics |2^, |3(j , and for studying the universal fluctuations 
in the impedance [j| and scattering matrix [3l| for a 
wave chaotic system. The cavity is driven by two ports; 
each of which consists of the center conductor (diame- 
ter 2a= 1.27mm) of a coaxial cable that extends from the 
top lid of the cavity and makes contact with the bottom 
plate of the cavity (Fig. 3(b)). From direct S21 measure- 
ments, we estimate that the cavity has a typical loaded 
Q of about 150 between 4-5GHz and about 300 between 
ll-12GHz. This translates to a damping parameter of 
A > 0.5 for the entire frequency range of this experiment. 
Hence we examine experimentally the time-reversal sym- 
metric (GOE) cases for the Z and ^-variance ratios in 
the high damping limit. 

To perform an approximation to ensemble averaging, 
two perturbers (shown gray in Fig. 3(a)), made up of 
rectangular ferromagnetic solids wrapped in Al foil (of 
dimensions 26.7 x 40.6 x 7.87 mm 3 ), are systematically 
scanned and rotated throughout the volume of the cavity 
by means of a magnet that is placed outside the cavity. 
The ensemble set consists of one-hundred different posi- 
tions and orientations of the perturbers within the cavity. 
The perturbers are on the scale of half the wavelength or 
bigger over the entire frequency range of the experiment. 
The full two by two S matrix is measured between 4 
and 12 GHz for each position of the perturbers. Below 
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FIG. 3: (a) The physical dimensions of quarter bow-tie 
chaotic microwave resonator are shown along with the po- 
sition of the two coupling ports. Two metallic perturbers are 
shown in gray, (b) The details of the coupling ports (anten- 
nas) and cavity height h are shown in cross section. 



4 GHz the mode density is too low to obtain meaning- 
ful ensemble averaging, while above 12 GHz the coupling 
becomes too weak to couple to the modes of the cav- 
ity, at least for this port geometry. Once the S matrix 
has been measured, it is then converted to Z through 
Z = Z (I + S)(I- s)- 1 . 

To eliminate the effect on the average of short ray 
orbits returning to the antenna (these lead to rapidly 
varying-with-frequency systematic deviations of the av- 
erage from the ensemble average, as discussed in p^ 1 ) 
we perform frequency averaging over a sliding window 
of width 300 MHz. We denote such sliding averages of 
impedance and scattering variance ratios by VR Z and 
VR S . The inset in Fig. 4 shows VR Z (solid line) over a 
frequency range 4 GHz to 12 GHz. Denoting the average 

of VR Z over the entire range, 4GHz-12GHz, by VR Z , we 
obtain VR Z = 0.49, and we find that \VR Z -VR Z \ < 0.02 
over the entire frequency range. This value of experimen- 
tally obtained VR Z is close to the ideal theoretical value 
of 1/2 for large damping. Also shown in the inset is 
the variance ratio obtained with no frequency averaging 
(small circles) These are deviations from the frequency 
averaged ratio values with a standard deviation of 0.04. 
Nevertheless the mean value of the variance ratio over 
the entire frequency range is 0.49. 

The circles, stars and dashe s in Fig. 4 show the 
variation in In [Var [Z21]], hi[y/Var[Zii]Var[Z22]] and 

respectively, as a function of 
the dashes 



lnfFi^ y/VarlZulVarlZn] 
frequency. The agreement is quite good (i.e. 
overlie the open circles) at all frequencies. 

Similarly in Fig. 5, we present data for the scatter- 
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Frequency (GHz) 



FIG. 4: Var[Z 2 i] (circles), ^"Var[Zu]Var[Z 2 2](sta,rs) and 

VR Z \J V^arJZn] Var [Z22] (dashes) are plotted on a natural- 
logrithmic scale as a function of frequency from 4 to 12 GHz. 
Inset shows the ratio VR Z versus frequency (solid line). The 
small circles show VR Z without any frequency averaging. 

ing variance ratio. Experimentally we obtain VR S — 
0.50 and \VR S - WW S \ < 0.08 over the range 4 GHz 
to 12 GHz. The circles, stars and dashes in Fig. 5 
show the variation in Var [S21], \/Var[Sii]Var[S22] and 
VR S \J Var [Si 1 ] Var [S22], respectively, as a function of 
frequency. Similar to the impedance data (Fig-EJ, we ob- 
serve that the data for V R S ^J Var [Sn]Var[S 22] (dashes) 
overlie the data for Var [621] (open circles), again indi- 
cating that the experimentally obtained value for VR S 
shows good agreement with the asymptotic theoretical 
values for highly damped time-reversal symmetric sys- 



tems over a large frequency range. 




Frequency (GHz) 



FIG. 5: Var [S21] (circles), ^/Var[Sii]Var[S 2 2] (stars) and 
VR S \/V ar[Su]V ar[S 2 2]{dsshes) are plotted as a function of 
frequency from 4 to 12 GHz. Inset shows the ratio VR 3 versus 
frequency. 

To sum up, we have used random matrix theory in con- 
junction with the radiation impedance characterizing the 
system-dependent coupling details to evaluate the vari- 
ance ratios (Q and JSJ. The main result is that the 
impedance variance ratio P is a universal function of 
the loss in the scatterer. 
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J. P. Parmantier who pointed out to us the result of 
Ref. [3 • We also acknowledge discussions with P. A. 
Mello and D. V. Savin. This work was supported in part 
by the DOD MURI for the study of microwave effects 
under AFOSR Grant F4962001 10374, and DURIP grants 
FA95500410295 and FA95500510240. 
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